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Magneto-transport of hard core bosons (HCB) is studied using an XX Z quantum spin model 
representation, appropriately gauged on the torus to allow for an external magnetic field. We find 
strong lattice effects near half filling. An effective quantum mechanical description of the vortex 
degrees of freedom is derived. Using semiclassical and numerical analysis we compute the vortex 
hopping energy iv, which at half filling is close to magnitude of the boson hopping energy. The 
critical quantum melting density of the vortex lattice is estimated at 6.5 x 10~^ vortices per unit 
cell. The Hall conductance is computed from the Chern numbers of the low energy eigenstates. At 
zero temperature, it reverses sign abruptly at half filling. At precisely half filling, all eigenstates are 
doubly degenerate for any odd number of fiux quanta. We prove the exact degeneracies on the torus 
by constructing an SU(2) algebra of point-group symmetries, associated with the center of vorticity. 
This result is interpreted as if each vortex carries an internal spin-half degree of freedom ('vspin'), 
which can manifest itself as a charge density modulation in its core. Our findings suggest interesting 
experimental implications for vortex motion of cold atoms in optical lattices, and magnet-transport 
of short coherence length superconductors. 

PACS numbers: 05.30.Jp, 03.75.Lm, 66.35.-|-a, 67.85.d 



I. INTRODUCTION 

Hard core bosons (HCB) are often used to describe 
superfluids and superconductors which are characterized 
by low superfluid stiffness and short coherence lengths. 
As such, HCB are relevant to cold atomic gases in op- 
tical lattices [1, 2], low capacitance Josephson junction 
arrays [3-5], disordered superconducting films [6], and 
cuprate superconductors [7-13]. 

At low densities, HCB can be treated by weak cou- 
phng (Bogoliubov) perturbation theory [14]. Closer to 
half filling, lattice umklapp scattering and the hard core 
constraints become important. Recent calculations of 
the dynamical conductivity of HCB near half filling [15] 
demonstrate the breakdown of weak-scattering Drude- 
Boltzmann transport theory in this regime. HCB exhibit 
so-called 'bad metal' phenomenology, (i.e. large resistiv- 
ity, linearly increasing in temperature). Such behavior 
has been often observed in unconventional superconduc- 
tors [16]. 

This paper also concerns dynamical correlations of 
HCB and their vortices near half filling. These will be 
exposed by including a weak orbital magnetic field in the 
Hamiltonian and studying the Hall effect. 

Our primary results are as follows. Firstly, we apply a 
combination of semiclassical analysis and exact diagonal- 
ization to the gauged XX Z Hamiltonian on a finite latice 
on the torus. We highlight the (sometimes overlooked) 
fact that a uniform magnetic field of one flux quantum 
penetrating the surface of the torus beaks translational 
symmetry. As a consequence, the semiclassical vortex 
center is subjected to a confining potential minimized 
at a well defined position. Fitting the low many-body 
spectrum to an effective single-vortex Hamiltonian, we 



determine the vortex hopping rate (effective mass). 

Near half filling, the vortex mass is found to be similar 
in magnitude to the HCB mass. This allows us to esti- 
mate the critical field for quantum melting of the vortex 
solid (superfluid) phase at 6.5x10^* flux quanta per unit 
cell. Secondly, at half filling we find doublet degeneracies 
associated with an odd number of magnetic flux quanta 
penetrating the torus. We associate them with symme- 
tries about the vortex position, and label the emergent 
degrees of freedom as 'vortex spin' (v-spin). Physically, 
these degrees of freedom correspond to the orientation of 
the charge density wave in the vortex cores. 

Finally, we compute the Hall conductivity using ther- 
mally averaged Chern numbers. In stark contrast to 
continuum bosons, and to electrons in metallic bands, 
we find that the Hall conductivity of HCB reverses sign 
abruptly at half filling. The associated Hall temperature 
scale vanishes at half filling, signaling a possible quantum 
phase transition for the thermodynamic system in a mag- 
netic field. Some of these results were briefly reported in 
a recent Letter [17]. 

This paper is organized as follows. In Section II the 
HCB Hamiltonian is introduced, with a discussion of its 
charge conjugation symmetry about half filling. Semi- 
classical approximations are derived in Section HI, for the 
various regimes of filling. At low density, we recover the 
Gross-Pitaevskii theory with its Galilean invariant vor- 
tex dynamics and classical Hall effect. At half filling, the 
continuum limit corresponds to the anisotropic gauged 
non linear sigma model. Its vortices possess localized 
charge density waves in their cores. Section IV describes 
the mathematical peculiarities of the gauged torus, in- 
cluding translational symmetry breaking (elaborated in 
Appendix A). Definitions of null lines, null points and 



vorticity centers are provided. The point group sym- 
metry generators Ily and ily are constructed and their 
commutator is calculated. The proof of v-spin degenera- 
cies at half filling is provided. Section V computes the 
vortex effective hamiltonian by combining semiclassical 
and exact diagonalization calculations. The critical field 
for quantum melting of the vortex lattice is deduced from 
our value of vortex hopping rate. Section VII computes 
the Hall conductance on the torus as a function of den- 
sity and temperature. We conclude in section VIII and 
discuss experimental implications of our results in cold 
atoms and cuprate superconductors. 



II. HARD CORE BOSONS 

The conventional Bosc Hubbard model for interacting 
bosons is 

(ij) i 

+iC/^n,(n,-l) , (1) 

In the hard core limit U — > oo, Eq. (1) reduces to the 
HCB Hamiltonian as "H = VHu^qV, where V is the 
projector onto the subspace where n^ = or 1 for each 
site. 

We use units where h ~ c ~ 1. (ij) denotes a nearest 
neighbor link on the square lattice; the lattice constant is 
a = 1. J is the Joscphson coupling, q is the boson charge 
and A^j the electromagnetic gauge field on a bond. T^ is a 
nearest neighbor repulsive interaction. In the HBC limit, 
The chemical potential fj, = corresponds to a density of 
half filling (n) = ^ , with half a boson per site on average. 

As is well-known, HCB operators obey an algebra cor- 
responding to spin-i; 



«1 = 


= Va\v = S+ 


a, = 


= Va^V^S- 


^i - 


= ~a% ^St + \ 



(2) 



By [Sf^Sj ] = 2Sf 5^p HCB operators obey constrained 
commutation relations. 



[a, , a]] = (1 - 2n,) 5, 



(3) 



The constraint effects of —2n^5ij become important near 
half filling, limj/^oo Hu is thus represented by the 
gauged spin-half quantum XX Z model, 

m 

+4l^^5f5|-ME(^^' + 5)- (4) 



It is widely believed that the ground state of Eq. (4) 
exhibits magnetic order. In the regime of T^ <C J, which 
is relevant to this paper, the ordered moment lies in the 
XY plane, {Sf) ^ 0. That is to say, except for the 
limits n = 0,1, the ground state of HCB exhibits long 
range superfluid order. 



A. HCB charge conjugation symmetry 

Another important distinction between the HCB 
Hamiltonian (4) and the finite U Bose-Hubbard model 
of Eq. (1), is the emergence of charge conjugation sym- 
metry in the infinite U limit. One defines the unitary 
charge conjugation operator, 



C = exp(z7rY,Sf 



(5) 



C transforms "particles" into "holes", i.e. C^h^C = 1 — h^, 
and 



C^n{qA,fi)C ^n{-qA,-fi) 



(6) 



At half filling (/i = 0), and A = 0, the Hamiltonian is 
invariant under charge conjugation on any lattice struc- 
ture [18]. 

A consequence of (6) is that the Hall conductivity 
(which is linear in q) is antisymmetric in the deviation 
from half filling, ie 



CTxyin, T) = -CTxyi^ - n, T) . 



(7) 



In contrast, the superfluid stiffness Ps{n) and longitudi- 
nal conductivity axxi^) are symmetric under n — >■ (1 — n). 
In terms of vortex motion, (7) implies that below and 
above half filling vortices drift in opposite directions rel- 
ative to the particle current. 



HI. SEMICLASSICAL THEORY 

The partition function of HCB can be represented by 
the spin half coherent state path integral [19, 20], 



Z = Vfi{T) exp / dr [iK - H' 
where 



(8) 



(9) 



/J oii-± vj v,wo (^Y'l v^j ' y^-^ij 



gA„) 



H''\fl,A\ = - J y^ sin t/; sm t/j cos ( 

-h y E cos Q, COS ^j - - X] ^^^ ^i- (1^) 



^i = {Oi, (j^i) are the polar angles on a sphere. The spin 
size S = ^ plays the role of the large parameter which 
controls the semiclassical expansion. 

In the classical (saddle point) approximation, for A = 
0, the ground state supcrfluid stiffness is 



d^H"" 



Ps = <1 



^\,r+x 



A=0 

= j(sin0r) = 4J7i(l - n). 



(11) 



which (in contrast to continuum bosons) exhibits a non- 
monotonic dependence on n. At half filling (optimal 
density), p^ is maximized. Quantum corrections en- 
hance pf{n = i) further by about 7% [21, 22]. The 
supcrfluid stiffness vanishes at the Berezinskii-Kosterlitz- 
Thoulcss (BKT) [23] transition temperature, computed 

to be Tbj^t - 1-41J. 

The kinetic term K of Eq. (8), determines the quan- 
tum dynamics. The harmonic spin- wave expansion of 
(10) yields a linearly dispersing phase fluctuations mode. 
The order parameter is suppressed to zero at all finite 
temperatures, in accordance with the Mermin- Wagner 
theorem. 



A. Low density, Gross-Pitaevskii limit 

For large negative values of the chemical potential ^, 
the action in Eq. (10) can be expanded around the fer- 
romagnetic (low density) state of 9i w tt. 



cosOi — > 2ni — 1 , sin^i « 2^ni{rn — 1) . (12) 
We define the continuous field 



V'(a?i) = 



,i<Pi 



(13) 



where a is the lattice constant, and replace the measure 
by 

Jl V cos Oi V(j>i — > X>^* PV n ® I 1 ^ A^^ ^*'^ ' 

(14) 
up to an unnecessary normalization constant. The Heav- 
iside functions enforce the hard core constraint n^ < 1 
in the i**^ unit cell at each time slice. In the low density 
limit, these constraints are ignored, and the action (10) is 
expanded to leading order in n^ ^ i and gradients \/ip. 
This yields an effective Gross-Pitaevskii (GP) theory [1], 



^GP 



5, 



GP 



y"pv*2?vcxp(-5Gp[V'*,v,^] 



(fx Idt 



+ 2^\^-'^~'^^^'^\ +25lV^l'. 



(15) 



where the effective mass and interaction parameters are 
given by 



1 

1* = 

16J 

g = 16(J + y). 



(16) 



In the presence of a magnetic field Bz, a density of 
TT-i, = B/(J)q vortices is produced, where 0q = ^ir/q is 
the flux quantum. The core profile function /(r — Rj) 
near vortex j is well approximated by minimizing the GP 
energy, which yields [1] 



fGp(r) 



ve 



(17) 



where ^ = l/^gm*n is the coherence length. For n ^ ^ 
one has ^ ^ a. The core density depletion is proportional 
to 1 — I/gpP- Hence it decays as 1/r^ away from the 
vortex center. 

In the high density limit, ri — > 1, the partition function 
can also be approximated by the same GP action (15) 
following a particle hole transformation (6). In this case, 
IV'P represents the density of holes. 

By neglecting the higher order gradients, and the hard 
core constraints, the GP theory does not include lattice 
scattering effects as it is completely Galilean invariant. 
Consider an externally induced uniform current density 



J = qnvs 



(18) 



In the moving frame of velocity Vs the vortices are sta- 
tionary. Therefore, back in the lab frame, a purely trans- 
verse electromotive field is produced by the moving vor- 
tices. 



-z X 3v 
—n^z X V-y 

q 

= 72— ^^■?- 



q^ n 



(19) 



That is to say, in the pure GP theory, the longitudinal 
(dissipative) conductivity vanishes, and the Hall conduc- 
tivity equals to the classical value. 







'xy 



nq 



(20) 



Spoiling Galilean invariance by the presence of nonuni- 
form potentials, boundary conditions, or by an underly- 
ing lattice can allow vortices to tunnel between different 
real space positions, resulting in a longitudinal conduc- 
tivity [24, 25]. 



B. Half filling, anisotropic cr-model 

Toward half filling, lattice scattering modifies the vor- 
tex structure and dynamics. At half filling /x = 0, 
the semiclassical theory of Eq. (10) is described by the 
anisotropic Non Linear cr-Model (NLSM) [26]. After a 
sublattice rotation OiGS e*'^'^* all the pseudo-spin inter- 
actions are antiferromagnetic. The spins fli are repre- 
sented by 



a, = rj, h{x,) ^1 - {L{x^)/Sy + L{x,)/S , (21) 

where rj^ = ±1 on the A {+) and B (— ) sublattices, 
respectively. The Neel vector h satisfies n^ ~ 1 and is 
orthogonal to the local magnetization L, i.e. h ■ L = 0. 
The complex combination 

n^=n,+*n,^ = |n^|e^^ (22) 

defines the local superfluid order parameter, and n^ cor- 
responds to a bipartite charge density wave (CDW) with 
two possible signs. Following Refs. [19, 20, 26], we sub- 
stitute Eq. (21) in the measure and action of (10), and 
expand them to quadratic order in L and Vri.. Integrat- 
ing out L arrives at the anisotropic NLSM path integral, 



'^NLSM — / ^ne e , 



(23) 




FIG. 1: Illustration of three static vortices with their v-spins. 
Arrows describe the v-spin directions, their ^-component is 
the local charge density wave near the vortex core, illustrated 
as ripples in the surface. Current density is depicted by black 
field lines. 



Vortex configurations at half filling can be viewed as a 
localized meron (half skyrmion) of the Neel field. Since 
\nxy\ = at the vortex center, and n1 = 1 — \nxy\'^, the 
semiclassical vortex has a CDW in its core, as illustrated 
in Fig. 1. 

Due to the finite anisotropy 'mass' m^ > 0, n^(r) de- 
cays exponentially away from the center 



nz{r) 
6 



v/l-n2(r)e'^(-), 



(28) 



where S-^ = Jdrfd-x £g is the Euclidean action, with 



^E = iX±|«_l| +5X>z 



(24) 



-^pA{V-tqA)nA^ + yi {VnJ^ + mlnl 



and 



T = 5'^?7i /dr (l-n^). 
^ 



(25) 



is the contribution from the geometric phase. The bare 
coupling constants are obtained directly from H^^: 



^^ 8J 



4.iJ + V)' 



and 



J 



V 



2{J~V) 



(26) 



(27) 



For A ~ 0, the isotropic (Hcisenberg) limit is at J = 
V, vriz — 0. The Neel ground state implies degener- 
acy between superfluid and CDW order, and the exis- 
tence of two massless Goldstone modes. At finite XY 
anisotropy, niz > 0, and there is one massless (phase) 
mode, and a gapped CDW (roton) mode at the CDW 
ordering wavevector (tTjTt). 



Indeed variational calculations have previously shown 
that at half filling CDW ordering is found in the local- 
ized vortex core [27]. In Section IV we shall show that 
the 'orientation' of the charge density wave is actually a 
continuous SU(2) symmetry of the quantum Hamiltonian 
at half filling, which we name v-spin. 

Since the system is charge conjugation symmetric at 
half-filling, there is no net charge depletion associated 
with the vortex core, and thus the statistical Berry phase 
for exchanging two vortices is zero. In other words, the 
vortices exhibit mutual Bose statistics. This is to be con- 
trasted with GP vortices at low filling. As shown in (17), 
GP vortices involve a large density depletion (or accumu- 
lation, above half filling), which decays slowly away from 
their core [28]. 

In the limit where the number of lattice sites TV tends 
to infinity, the confining potential on the vortex vanishes, 
and the vortex energy is periodic on the lattice. Its min- 
ima lie in plaquette centers {i.e. at dual lattice sites). 

When a vortex moves between dual lattice sites, the 
path dependent geometric phase T yields 27r times the 
number of bosons enclosed by the path. At half filling, 
this amounts to an effective tt flux per dual plaquette. 
These phases can be incorporated in an effective hopping 
model by the dual lattice gauge field Aj^ j^, along the 
link from site i? to R + rj. Thus for a single vortex 
on the infinite lattice, one can write an effective Harper 



Hamiltonian, 

H^ = -¥vY. i^'^'"'''' ^R^R+v + H-^- 
R,v 

E -^RM+v = 2^^ E "- (29) 

C ieint(C) 

where the sum on the second hne is a over a set of Unks 
comprising a closed path C on the dual lattice, and int(C) 
is the interior of this path, which consists of a set of sites 
on the original lattice bounded by C. 



IV. THE GAUGED TORUS 

We now return to the original HCB Hamiltonian, 
Eq.(4). We consider a finite square lattice, of dimensions 
L^ X Ly, with N = L^Ly sites and periodic boundary 
conditions in both the x and y directions. This toroidal 
geometry is convenient for the study of finite lattices as 
it minimizes the effects of boundaries. It also provides 
external control over the positions of vortices via the two 
Aharonov-Bohm (AB) fluxes which run along the two cy- 
cles of the torus. The lattice site positions are labelled 
as 



Xi 



0,1,2,, 
0,1,2,, 



■ 1 J^x 

■ 1 J-^u 



(30) 



A uniform magnetic field B is everywhere perpendicular 
to the surface, such that the total number of flux quanta 
penetrating the surface is Nj = NB/cpQ, where (J)q = 
2'!T/q is the flux quantum. 

We construct a (piecewise differentiable) gauge fleld 
A{x) which interpolates the lattice gauge field on the 
surface of the torus, and obeys 



z-V X A 
r]-A{r + ^T]) 



B 



T+^ ' 



(31) 



where 77 = ±x, ±y. A determines the magnetic fluxes 
which flow through vertical and horizontal circumfer- 
ences of the torus. These are given by the gauge invariant 
Wilson loop functions, 

Wy{x) = q i>dy Ay{x,y) mod 27r 

W^{y) = q (f)dxA^{x,y) mod 27r. (32) 



The dimensionless AB parameters = {0^,0 ) are de- 
flned by the Wilson loops at a; = and y = 0, 



0y = Wy{X = 0) 

O. = W^y^Q). 



(33) 



lives on the reciprocal torus [0, 27r) x [0, 27r). 

W^ = 0,a = x,y define the null lines on the torus. For 
Nj, = 1, there is one null line in each direction x = Xq, 




FIG. 2: The gauged torus, which defines the finite size geome- 
try of the HCB. The torus surface is penetrated by a uniform 
magnetic field of one total fiux quantum, and threaded by 
two Aharonov Bohm fiuxes = (Q'^jQ^). Thick (red on- 
line) circles denote the null lines which enclose zero flux, and 
intersect at the null point Ro{0). The vorticity center i?v 
is located on the antipodal point to Ro- The circulating cur- 
rents of the ground state are illustrated by thick green arrows. 
This geometry is used to compute the vortex mass and Hall 
conductivity of HCB, and to prove the v-spin degeneracies at 
half filling. 



and y = Yq, as depicted in Fig. 2. Their intersection is 
the null point Rq = (Xq,Yq), which constitutes a gauge 
invariant symmetry point on the torus. 

Xoi&) = - 



Yo{0) = 



2tt 

Ly 0,J. 

2^ 



(34) 



The existence of a special point Rq on the torus, demon- 
strates the unintuitive fact that a uniform magnetic 
field necessarily destroys lattice translational symmetry. 
This fact is closely related to the quantization of Dirac 
monopoles in three dimensions. We elaborate further on 
this fact in Appendix A. Eq. (34) shows that Ro{0) can 
be moved continuosly on the torus by changing the AB 
parameters [29, 30]. 

As we shall see in Sec. V, semiclassical analysis and 
exact diagonalizations find that the center of vorticity 
Ry is located at the antipodal position of the null point 
on the torus, 



i?v(0) = (ii. 



Xom > ih 



Yoi0)) • (35) 



For larger magnetic fields, Ncf, > 1. there are A^^ null 
lines in each of the x and y directions. This introduces a 
set of TV? null points which form an evenly spaced square 
lattice (which may or may not coincide with the original 
lattice sites). These are indexed by m, n = 

1 



,N^-l 



Rr{0) = Roi0) + -^{mL, , nLy), 



(36) 



(e,,e,) = (o,o) 



(e„,e,j)=(3ir/£,,37r/r.,,) 
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FIG. 3: The lattice gauge field on links, A.r,r+ri, according 
to Eq.(39), for two choices of AB parameters {Ox, Qy)- The 
length and thickness of the arrows are proportional to the 
magnitude of A. For a single flux quantum, A^^ — 1, the two 
null lines are marked by red lines. The null points Ro{0) 
and vorticity centers R\j{0) are depicted. Note that in the 
right flgure, both points are located at plaquette centers. 




FIG. 4: Emergence of V-spin' degeneracies. Four lowest 
eigenenergies of HCB on the torus, with N^ — 1, as a func- 
tion of the AB parameter 0'', with Q^ — 0. The spectrum 
separates into doublets. Notice the exact degeneracies which 
occur when the vorticity center coincides with a lattice posi- 
tion. 



Correspondingly there are TV? vorticity centers, 

Rr{0) = Rri0) + -2iL.,Ly). (37) 

A. Choosing a gauge 



with flux B ~ N^<J)q/ LxLy in each plaquette, and the 
two Wilson loop functions 



Wy{x) ^ xqBLy + Qy, 
W,{y) = -yqBL, + Q, 



(41) 



The uniform magnetic field must integrate to an inte- 
ger number of flux quanta N^, 



B 



2nN^ 



qLxLy 
The gauge field is given by 



(38) 



AX 

r,r+x 

AV 
r,r+y 



- mod{y-Yo,Ly)BL,H{Xo,x) 
mod{x- Xo,Lx)B . (39) 



Note that for < m < n, mod (— m, n) = n—rn. The 
function H{Xq,x) ensures that A^ ^ ^ vanishes unless r 
is immediately to the left {—x) of the null line x = Xq. 
It is defined by 



H{Xo,x) = 



1 0< mod (Xo-a;,^^;) < 1 

otherwise. 



(40) 

For a continuous position r we define A"(r) to be the lin- 
early interpolated gauge field between the two enclosing 
links in the a direction. 

For = (0, 0) the nuU point R„ is at R^ = (0, 0), and 
the vorticity center is therefore at R^ = {^Lx,^Ly). 

Our gauge choice is shown in Fig. 3. The gauge invari- 
ant content of A consists of the uniform magnetic field 



B. V-spin degeneracies 

In the process of calculating the Hall conductance (see 
Section VII), we computed the spectrum at half filling, 
for an even number of sites, with one total fiux quan- 
tum of magnetic field. We encountered a sequence of 
AB fiuxes 0i, where the whole spectrum becomes two- 
fold degenerate. These degeneracy points are demon- 
strated in Fig. 4 for N^ = 1, for the lowest two mul- 
tiplets. The level crossings indicate the existence of a 
non-commuting symmetry generators [31], which act on 
the wavcfunctions of vortices introduced by the external 
magnetic field. We now construct these symmetry oper- 
ators and compute their commutation relations. 

As discussed earlier, for a finite magnetic field (TV^ > 
0), Ji does not possess the lattice translational symmetry. 
Nevertheless, with respect to the vorticity center Ry{0) 
we can define two reflection operators 



P^{x,y) 



(mod (2Xv - 
(x, mod {2Yy 



■,Lx),y), 
y,Ly)), 



(42) 



which by (35) are equivalent to reflections about Rq. 

Now, by appropriately tuning using Eq. (34), the 
vorticity center Ry can be chosen to coincide with a sym- 
metry point of the square lattice, such as any lattice site. 
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FIG. 5: Figuring out Nnuii for two cases of even size latices. 
The parity of the number of lattice sites (black dots) which lie 
precisely on the null lines (red color online), determines the 
commutation rule of 17^ and U^, according to Eq. (53). The 
two examples explain why N„uii is odd whenever the vorticity 
center i?v is positioned precisely on a lattice site. 



bond center or plaquette center. Reflecting the Hamilto- 
nian about that symmetry point, leads to 



Xct A 



(43) 



The gauge invariant content of A" describes an inverted 
uniform magnetic field B = —B, and a reversed sign of 
the Wilson loop functions (32). 

The reversal of the fields in A can be undone, at half 
filling, by applying the charge conjugation transforma- 
tion C (5), and a pure gauge transformation U". Thus, 
we construct two operators. 



where. 



and 



ny = uycp^ 



i7« = exp(z^x"(r)5|i) , 



x"(r) = ydr'.(A(r') + A"(r')). 



(45) 



(46) 



In the line integral we use the interpolated gauge field 
defined after Eq. (39). Since A and —A" describe the 
same magnetic fields, they obey. 



V X L4 + A" =0 



(47) 



This implies that x" is independent of which continuous 
path (of zero winding number) is chosen between Rq and 



It is easy to verify by this construction that for all 0i 
such that Ry{@j) is a symmetry point the 77y operators 
become symmetries of the Hamiltonian: 



[H[A],i7^]=0 , a^x.y. (48) 

Now we calculate the commutation relation between 
TTy and Ily . This is a straightforward but slightly tedious 
procedure. Using the gauge choice (39), and (46), 



X'ir) = ^^mod(y-yo,iy)(l-'5^,Xo) 

Lly 

X'{r) - 0, 



(49) 



Note that Xy vanishes on the null lines. Multiplying the 
two LLy operators yields 



r 

n<^ ny = exp U ^ (x^ - x^ (^v" [r]))s^) py P^ 



= e-^^7T^7T^, (50) 

where we have used [Py , Py] =0. The overall phase is 
given by the operator 

r 

oor = x^ + x"(^v'M) -x"- X^(^v M) . (51) 

It can be directly verified from (49) that 



r G null lines 

2ttNa, otherwise. 



(52) 



(•44-) Since exp (i2TT'mS^) = {—lY 



g-.T ^ (_i)iV,(JV-iV„„n) ^ ^_if,N„^u^ (53) 

where A^nuii is the number of sites which sit precisely on 
the two null lines. For even N^, e*^ = 1 and [iTy , ily] — 
0. 

For odd N^f,, and odd A^nuii, one obtains e~*^ = —1. 
Let us prove a simple lemma concerning the parity of 

iVnull. 

Lemma: For even size lattices, if Ry is tuned to be 
precisely on a lattice site, then Nnuii is odd. The proof 
is illustrated in Fig. 5. 

Proof: Since we assume that N is even (to describe 
precise half filling), there are two cases to consider: 
(i) For an even by even lattice, L = (2m, 2n), if we choose 
Rq on a lattice site, it is easy top see that Ry must also sit 
on a lattice site. The number of sites which contribute to 
Nnuii arc the sum of lattice sites in the x and y directions 
minus the null point itself which is counted twice: 



iV^;;ii = 2m + 2n-l. 



(54) 



Hence e'^ = — 1. 

(ii) For the odd by even lattice, e.g. L = {2m + 1, 2n), 
we choose i?Q, to be in the middle of a bond in the x 
direction. The null line includes only the sites on the 
x-nuU line which is odd: 



ATCO 



2m + 1 



(55) 



Thus, here too e~*^ = -1. Note that in both (54,55), the 
vorticity center is situated on lattice sites R^ = {m,n). 
QED. 

Thus we conclude that for an odd number of fluxes 
Nij, = 2n+l,ii Ry{0i) is located precisely on any lattice 
site, then Uy and 11^ anticommute. 

Under these conditions, all states of T-L[0i\ must be 
at least two- fold degenerate. This follows the standard 
proof: Since 



[H,7T^] = 0, 



(56) 



and Uy has eigenvalues ±1, then each common eigenstate 
of % and iJy, can be labelled by |i?„, tt^ = ±1). Now, 



iT^iT^|£;„,l) 



-i7^iT^|S„,l)«|£;„,-l), (57) 



that is to say each eigenenergy En is associated with 
a degenerate pair of eigenstates with opposite quantum 
numbers of Uy. 

Uy and Ily are point group symmetries about the vor- 
ticity center. We can also construct a third symmetry 
operator H^ as 



m 



-in^nl 



(58) 



The three operators IIQ = {Uy , Ily , Ily) are unitary 
and Hermitian, 

ne = (ne)t = (ne)-i = 



[n^Y = 1. 



(59) 



Therefore their eigenvalues are ±1. The II operators be- 
have as Pauli matrices and can be used to construct an 
SU(2) algebra of spin half. 



T-" — in" 



a = x,y,z 



(60) 



We note that for multiple number of magnetic fluxes 
N^ > 1, degeneracies appear for any odd number of vor- 
tices. This is consistent with the Kramers doublets asso- 
ciated with an odd number of interacting spin half par- 
ticles. 



where n were defined in (21) and Da are gauge invariant 
derivatives in the xy plane. In the absence of a gauge field 
A = 0, a single meron (half a skyrmion) of a continues 
classical field has topological charge 



^-h 



dVT^(r,A = 0) = ±i . 



(62) 



A HCB representation of r^ is constructed spin-half 
operators. In the presence of one flux quantum, wc define 
the modified topological charge operator as 



iT^rriA) 



(63) 



Q is not expected to be quantized at ±i. However, 
we have found that in the low lying eigenstates of (4), its 
sign correlates with the eigenvalues of H^ : 



sign((£;„,7r^|Q|£:„,7r^ 



(64) 



We conclude that tt^ of a single vortex measures the sign 
of the CDW in its core, with respect to sublattice A. 



V. THE HAMILTONIAN OF QUANTUM 
VORTICES 

A. Vortex confining potential 

The current density operator is given in the pseudo- 
spin representation by 



Jnir) = -2iJq {e''i^^'-+^S+S-+^ - H.c. 



(65) 



By choosing A to describe one flux quantum of uni- 
form magnetic flcld through the whole lattice one vortex 
is introduced into the low energy eigenstates. Indeed, 
we verified that the exact ground state exhibits a vortex 
pattern of the current density (j,j), defined in (65). Also, 
the center of vorticity is agrees with the value of Ry{0) 
as defined in Eqs. (34,35). 

We determine the effective confining potential on the 
vortex variationally. We choose the square geometry 
i^ = N, and define a vortex coherent state centered at 
i? by a spin coherent state |/7v[i?]). All the unit vectors 
fiy lie in the xy plane with the azimuthal angles (f>r given 
by a Jacobi theta function [32] : 



C. V-spin and meron density 

The semiclassical analysis of HCB vortices at half fill- 
ing, shown in Eq. (28) finds that the vortex has a CDW 
in its core. This is signalled by the local order parameter 
{n^{r)) 7^ 0, as illustrated in Fig. 1. 

We define the modified 'meron density' operator as 



T^{r) ~ h- Dxhy. Dyfi, 



Dr, 



Oct ^9-^a; 



(61) 



■r{R) = -Imlogi?(z(z-Z)-i-§ 
+7r(y-(i-y).x) 

d{z) = Y^ e-'^^'e^ 



-7rn Jli-nnz 



n— — oo 



(66) 

(67) 



Here we use scaled complex coordinates z — {x + iy)/ L 
and Z = {X -\- iY)/ L for the position and vortex center 
r and i?, respectively. 
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FIG. 6: A variational vortex configuration with its center 
located at i? = (0.5, 0.6)iy. The black dots mark the center 
of vorticity Jlv, where the variational energy is minimized, 
a) Phase field marked by directions of the arrows, b) Current 
density distribution where the colors correspond to the local 
current magnitude. Note the excess currents flowing around 
the torus in the negative x direction: a consequence of the 
vortex being displaced from variational minimum. 
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FIG. 8: Fitting the single vortex energies of the effective 
Harper Hamiltonian (70), to the many-body spectrum of HCB 
(4). The HCB Hamiltonian is defined on a 4 x 4 lattice at 
half filling, which is embedded on a torus which is penetrated 
by one fiux quantum. The parameters J = 1, and V = 
are chosen for this figure. Up and down arrows denote the 
v-spin magnetic quantum number r^. The confining poten- 
tial of the vortex in the Harper Hamiltonian is calculated by 
a variational calculation, Eq. (69). The effective vortex hop- 
ping rate is fit to be tv = 4J. The lowest three doublets of 
the two models agree within 2% of the first energy gap. 



FIG. 7: The energy as a function of vortex position, for a 
single vortex on the torus of dimensions L x L. The vorticity 
center is located at i?v 



\L. 



While <f>r is discontinuos on the torus at x,y = L, the 
gauge invariant current density 



{il^\jri\f2y) = gJsin 



\-r] -^r.r+rj) 



(68) 



We can now combine the single vortex hopping terms 
of Eq. (29) with the confining potential to obtain the 
Harper hamiltonian on the finite torus 



ttN 



2^y 



E 

R,r) 



'^'^^■^+''&W+, 



~UN{R)b^jibji. 



H.c. 



(70) 



is continuous. An example for the phase and current 
distributions is given in Fig. V A. 

The effective confining potential on the vortex is given 
by the classical energy defined in (10). By semiclassical 
estimates [33], its curvature at R^ scales as 1/A^. For 
lattices of size L > 4 we fit the variational potential by a 
two dimensional quadratic function, which scales as 1/N 
as, 

Un = ^K\R - Ry\'^/N, K = 39.2 J. (69) 

U is minimized at the vorticity center iiy, which was 
defined in (35). Fig. 7 depicts the confining potential as 
a function of vortex center for the choice of = (0, 0). 



B. Vortex hopping Amplitude 



For a quantitative quantum theory of vortices we need 
to evaluate the effective hopping ty. Since vortex tunnel- 
ing between lattice sites depends on short range many- 
body correlations, we extract ty from exact numerical 
diagonalizations of "H on 16 — 20 sites clusters, in the 
presence of a single flux quantum. 

By tuning t\j, we fit the lowest three eigenenergies En 
of T-L to those of the effective Harper Hamiltonian (70). 
The fit is shown in Fig. 8. 

Our primary concern is that the low eigenstates will 
not be exclusively described by H^ . since there are also 
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a) HCB vorticity density 
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b) Vortex wavefunctions probability density 
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FIG. 9: (a) The HCB vorticy density pi"\R) as defined in 
Eq. 72), for the states |n,t) whose spectrum is depicted in 
Fig. 8. A uniform background vorticity has been subtracted, 
(b) ProbabiUty density of the lowest single particle states of 
the Harper model given by Eq. (70). The qualitative similar- 
ity between (a) and (b) supports the fit of the vortex hopping 
rate, and the validity of the Harper Hamiltonian for single 
vortex dynamics. 



low energy superfluid phonons (phase fluctuations) [34- 
37]. However, we can estimate the phonons lowest exci- 
tation to be gapped by the finite lattice with the energy 
scale 2ttJ/L, which is larger than the energies we have 
fitted to (70). 

Our results for ty{n^^,V/J), for A^ = 20 can be de- 
scribed by the fitting formulas, 



iv(n,0) 



AJ - 50.4 J [n - 
4J + 6F+10.8 



i) -f 5056J(7i-i) , 
J ' 



(71) 



The system parameters were varied throughout the range 
\n-\\< 0.2, and V/J < 0.5. We find that at half filling, 
the vortex hopping rate ty varies very little between the 
TV = 16 and A^ = 20 lattices. This indicates that the 
bare vortex kinetic energy is determined by short range 
correlations, and thus does not require a large lattice to 
be computed with acceptable accuracy. 

To further test the association of the many body lowest 
eigenstates with vortex center fiuctuations, we measure 
the vorticity density defined on the dual lattice, 



o(") 



{R) 



plaq.i^ 

E 

r,T} 



(*«|>,r,|*n) V, 



(72) 



for each of the low lying states |^n). We compare pv to 
the single particle probability density of the correspond- 
ing wavefunctions of H^ . As shown in Fig. 9, for the fit- 
ted value of ty, the corresponding distributions increase 
in width in a qualitatively similar fashion. This demon- 



strates that the low lying eigenstates of H correspond to 
quantum fiuctuations of the vortex position. 



VI. QUANTUM MELTING OF THE VORTEX 
LATTICE 

At half fiUing, for A^^ > 1 semiclassical evaluation of 
vortex interactions at distances larger than the core ra- 
dius is 



U, 



int 



~TrJlog{\R,-Rj\). 



(73) 



Integrating out the phonon fluctuations, produces an in- 
stantaneous logarithmic (2D Coloumb) interaction be- 
tween vortices, plus retarded (frequency-dependent) in- 
teractions [35, 36]. Since we are interested in the short 
wavelength fluctuations which are responsible for quan- 
tum melting of the vortex lattice, we ignore these retar- 
dation effects. 

Thus in the large lattice limit, at a flnitc magnetic field 
B, the multi- vortex quantum hamiltonian is given by the 
Harper boson plasma (HBP) 



(74) 



-TT J 2J rij n^ \og{\Ri - Rj 



At half filling the continuum limit of (75) can be taken 
as follows. The single vortex dispersion exhibits a two- 
fold degeneracy of the ground states of the Harper Hamil- 
tonian at TT flux per plaquette. This implies the degen- 



eracy Ef^ = E^ 



fc+(0,7r): 



at low |fe| << TT. Since we wish 



to expand the hamiltonian at long wavelength, we retain 
the degeneracy by the v-spin label s ^tii- The vortex 
effective mass is deflned as 



Mv = 






= iAv 



(75) 



This leads to the continuum spin-half Coulomb Bosons 
(CB) Hamiltonian for the vortices as half flUing: 



ncB = E 



i,s=ti 



2^/. 



7rJ^log(|R,-i2,|) 



iy^J 



B 



l-.v5:iR,i 



$0 



(76) 



At low vortex densities, interactions clearly must dom- 
inate over the kinetic energy, and vortices form a vortex 
lattice. This a superfluid phase with v-spin correlations. 

At a flnitc temperature T™, if we ignore the quan- 
tum effects of the kinetic term, the classical melting tem- 
perature T„i is independent of vortex density (magnetic 
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field) [38, 39]. (Any change in B can be absorbed by 
scaling Ri appropriately, leaving the classical energy in- 
variant). 

However, quantum fluctuations of the Colounib plasma 
Eq. (76), increase with the vortex density i3/$o, until 
quantum melting is reached at a critical vortex density 
Bci- This melting is analogous to that of spinless CB, 
studied by Magro and Ceperley (MC) [40] by diffusion 
Monte-Carlo. 

MC used the dimensionless parameter to describe the 
CB density n^, 



Tm^a^, 



(77) 



where Oq is the Bohr radius. We set the Bohr radius to 
be 



Oo 



■kJAL 



1/2 



(78) 



to match between the model of MC and our "Hqb of 
Eq. (76). 

MC found that below Ts « 12 the boson lattice un- 
dergoes quantum melting, i.e. they found the critical 
quantum melting density of 



1 



1447ra2 ' 



(79) 



Above this density, the CB looses translationally symme- 
try breaking and becomes a quantum liquid, which will 
be discussed in Section VIII. Using our values of My 
from Eqs. (71,75), this translates into a critical vortex 
number per lattice site of of 



< 6.5-7.9- 



X 10 '^ vortices per site. (80) 



This is a suprisingly low vortex density, above which a 
vortex quantum liquid (QVL) is created. 



VII. HALL CONDUCTIVITY 

We have shown earlier by Eqs. (20 , 7) that the Hall 
conductivity in the low and high density obey the effec- 
tive Galilean invariant limits. 



<»^ 




FIG. 10: Hall conductance as a function of boson number TVf, 
for hard core bosons on a 4 x 4 lattice on the torus with one 
penetrating flux quantum. Temperatures vary in intervals of 
AT = 0.05J. The jump of the zero temperature conductance 
at half filling, is smoothened at finite temperatures. 



The zero temperature Hall conductance of a finite lat- 
tice embedded on a torus is defined by the Chern num- 
ber [41]: 



2-K 27r 

...(^)^|^/de./.e,im(||£ 



a^o 



99, 



(82) 



where l*o(0)) is the exact ground state of (4), in the 
presence of AB fluxes @. q^/h is the quantum of con- 
ductance. In the absence of degeneracies and level cross- 
ings, the Chern number axy{N)h/q^ is an integer. We 
compute Eq. (82) for a sequence of finite lattices. In Fig- 
ure 10 we plot (JxyiN) for a square lattice of size 4x4 
with aa single fiux quantum N^ = 1, as function of boson 
numbers Ni, = [0, 1, ... , 16]. We find that the Hall con- 
ductance follows two straight lines as given by Eq. (81), 
with an abrupt jump to zero at half filling. The same 
behavior was found for all smaller lattices, and refiects a 
sharp change in vortex dynamics around half filling. 

We extend Eq. (82) to finite temperatures by thermally 
averaging over all eigenstates \^n), 



nq 




n< i 


(l-n)g 
B 


1- 


-n< i 



(81) 



In terms of vortex motion, this relation implies that be- 
low and above half filling vortices drift in opposite di- 
rections relative to the particle current. In the following 
we shall study the transition between these two regimes. 
Since the continuum approximation is expected to fail 
near half filling, we resort to a numerical computation of 



CTxyiT) 



hn 



E h^Jde 



n=0 



xlm 



'E„/T 



d^n 



ae. 



d^n 



de. 



(83) 



En{0) and |^„(0)) are the exact spectrum and eigen- 
states of (4). The results obtained with Eq. (83) are 
matched at high temperatures with the the conductiv- 
ity calculated using the Lehmann representation of the 
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FIG. 11: axy of hard core bosons, as a function of temperature 
(in units of iv), for 4 to 7 bosons on a 4 x 4 lattice. For 
low temperatures, we calculate Oxy using Eq. (83), while for 
T > 1 we use Eq. (84). The dashed line is an interpolation 
between the two calculations. The points where a^y drops to 
half its value at T = are indicated. The inset shows the 
temperature scale Th as a function of density difference from 
half filling. 



Kubo formula [42] 



e-/3B„ _ g-/3B„ 



iO+ NZUJ ^-^ En - Em + LO 



lim 



■E 



x(vl/„|^j.(r)|vl/„)(vl/,„|^j,(r')|*„), 

r r' 

(84) 

where the current operator jr^ is defined by (65). 

The Kubo expression for axy is evaluated at high 
enough temperatures where the w — > limit is well be- 
haved. In Figure 11 we plot axy as a function of tempera- 
ture at different HCB densities, by interpolating between 
Eq. (83) at low temperatures and Eq. (84) at high tem- 
peratures. 

We see that in general, the magnitude of axy{T,ny^) 
decreases with temperature, and the discontinuity as a 
function of filling at zero temperature smoothens at fi- 
nite temperatures. As the temperature is lowered, the 
reversal of the Hall conductance takes place in a nar- 
rower region around half filling. A characteristic Hall 
temperature THin^) can be defined by. 



r^xyiTn) = ^CFxy{0)- 



(85) 



In the inset of Fig. 10, we show that Th increases with 
I'^b ~ 5I' although we cannot estimate the critical expo- 
nent from the small cluster calculation. 
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FIG. 12: Different vortex drift directions (purple arrows) in 
the presence of a bias current (red arrows), for regions of 
boson density n which is lower (blue online) and higher (green 
online) than half filling. 



VIII. DISCUSSION AND EXPERIMENTS 

In this paper, we have determined the vortex effective 
hopping rate iy, in Eq. (71), and the Hall conductivity 
axy, in Fig. 10 from the eigenstates of 16 site clusters 
on the torus. We emphasize that these quantities serve 
only as short wavelength "bare parameters", to be used 
to in the single vortex Harper Hamiltonian Eq. (75) and 
the multi- vortex Bose Coulomb liquid model (75). The 
charge transport coefficients depend on the thermody- 
mamic phases of the the latter model. 

Vortex Solid ~ The vortex solid phase, expected for 
vortex densities lower than n'^^ of Eq. (80), has superl- 
fuid (i.e. superconducting) order. Vortices are pinned 
to their lattice positions and therefore axx = 00 and 
(^xy = 0. The role of fy is to produce quantum zero 
point motion and effective v-spin super-exchange interac- 
tions, which are ferromagnetic (since vortices have Bose 
statistics). Below the v-spin ordering temperature Tvspin, 
charge density waves might be expected with significant 
magnitude in the vortex cores. This phase is a weak 
supersolid. Incidentally, density modulations have been 
observed near vortex cores of High Tc cuprate supercon- 
ductors [43] and analyzed within bosonic models [44, 45]. 

At low vortex density, the supercxchange interactions 
decay rapidly, which reduces Tvspin- For T > Tvspin the 
v-spins will contribute an extra entropy density given by 



'S'vspin = -— l0g2. 
<I>0 



(86) 



Quantum Vortex Liquid ~ At vortex densities which 
exceed 6.5 x 10"'^ vortices per lattice site, we expect 
the vortex lattice to melt and give way to the Boson 
Coulomb Liquid studied by MC [40] . Superfiuidity of the 
CB translates into a Mott insulating behavior of the orig- 
inal bosons [34]. However, MC have found that the liquid 
phase Eq.(76) is incomprcsible and hence exhibits vanish- 
ing condensate fraction [46]. Furthermore, retardation 
effects act to suppress dual superfluidity [47]. The value 
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of the transport coefficients of the QVL phase is there- 
fore left as an important open question. Away from half 
filling, our results for a^y show that the vortices are sub- 
ject to a strong magnetic field, which further suppresses 
their condensation. At low boson fillings and large vor- 
tex density, n^/n^ < 1, there is evidence for fractional 
quantum hall phases [48, 49]. 

The QVL phase discussed above is distinct from the 
vortex- antivortex condensate (VC) phases which were 
predicted at rational boson filling fractions, riy^ = p/q, 
in the absence of a magnetic field [45, 50-52]. These 
are expected at strong longer range interactions y « J, 
and corresponds to Mott-insulating commensurate CDW 
phases. 

Hall coefficient - The abrupt reversal of Hall coefficient 
was found for 16 and 20 site lattices. This effect corre- 
lates with the rapid change in the semiclassical vortex 
core profile at half filling, since we know by the GP equa- 
tion (17) that vortices have a diverging density depletion 
(accumulation) in the low and (high) filling regime, while 
they have a localized charge density wave in their small 
core at half filling. The sign of axy determines the drift 
direction of a vortex with respect to a bias current. This 
rapid reversal of axy may be relevant to the rapid change 
in Hall resistivity as a function of doping, which was ob- 
served in (Lai_:rSrj;)2Cu04 [53, 54]. 

In Fig. 12 we propose a set-up to observe the Mag- 
nus action reversal for cold bosonic atoms on a rotating 
optical trap [55]. If the density of bosons is allowed to 
vary slowly in space across half filling, we expect a rapid 
change of the vortex drift directions at the half filling 
line. The vortices would drift downstream with the bo- 
son current for n < ^ and upstream for n > i. 

Soft core interactions - We have not considered re- 
laxing the hard core constraints of the bosons, which 
would be described by the Bose Hubbard model (1). For 
U/J < oo, the charge conjugation operator ceases to be 
an exact symmetry at half filling, and Eq. (6) is not valid. 
Therefore the Hall coefficient will not be precisely zero at 
half filling, and the v-spin degeneracies will be lifted by 
the finite U/J corrections. A full determination of the 
Hall conductivity in the U/J versus n phase diagram, is 
an interesting open question. In particular, one would 
like to find out how the zeros of the Hall conductivity 
connect between the HCB limit and the free bosons limit. 
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FIG. 13: Stem's construction. A large number, A'^, of pos- 
itive magnetic monopoles (red circles) and N — 1 negative 
monopoles (blue squares) uniformizes the the radial magnetic 
field, as shown in Fig. 14. However, the internal fiux which is 
measured through the disks which cut through the torus, os- 
cillates wildly as a function of azimuthal direction x, as shown 
in Fig. 15. 



Appendix A: Translational symmetry breaking on 
the continuous torus 



In section IV we have shown that the ground state of 
the HCB Hamiltonian Eq. (4) with < N^ < N exhibits 
translational symmetry breaking (TSB) relative to the 
lattice periodicity, in both x and y directions. At first 
thought this is very surprising: the flux per plaquette is 
uniform, so one might expect that physical observables in 
any non-degenerate state would not distinguish one lat- 
tice point over another. However, one empirically finds, 
for example for N^ = 1, that the ground state current 
circulates around a preferred position, R^. 

In this appendix we first explain the reason for the 
TSB. By working out a specific proposal for trying to 
construct a purely radial magnetic field, we learn why 
TSB is in fact unavoidable on the torus. 

In the particular choice of gauge field in Section IV, 
the gauge invariant Wilson loop functions Wx , Wy de- 
fined in Eq (41) are linear functions of x and y modulo 
2tt. Therefore, by construction they break translational 
symmetry in both directions, and we use them to define 
the special null point on the torus, Rq. 

Now we show that by Stoke's theorem, if A is continu- 
ous on any interval [xi , X2\ , then Wy (x) must be picccwise 
linear, 



Wy{X2) - Wy{Xl) = qBLy{X2 " ^ 1 ) , 



(Al) 



where Ly is the circumference of the torus in the y di- 
rection. A{r) and Wy{x) however cannot be continuous 
everywhere on the torus, since Wy{x2) must be periodic 
for ^2 — > a;2 -l- Lx which is inconsistent with a continuous 
linear behavior given by (Al). 
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FIG. 14: The radial component of the magnetic field in the 
plane y = 0, containing the ring A*' = 20 monopoles and 
A'' = 19 anti-monopoles on the line y — 0, as shown in Fig. 13. 
as a function of the azimuthal coordinate x. The oscillations 
magnitude decays as a function of A'^ leading to a uniform 
radial field. 
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FIG. 15: Wilson loop function Wy{x), and the phase function 
C{x), for the monopole configuration of Fig. 13 for A'^ — 20. 
While Wy oscillates wildly at large A'' and does not converge 
to a limit, the physically relevant function ^ is continuous for 
all A'^ and breaks translational symmetry. 



The magnetic field of N^ = 1 enters the torus and must 
end in a magnetic charge. By Dirac quantization [56], 
'magnetic charge' density must be quantized as pm — 
9in 5{r - r,„) where. 



= 1/(29). 



(A2) 



Therefore the physical reason for the TSB is the point- 
like discreteness of the magnetic charge. Consider an 
embedding of the torus in three dimensions. Dirac quan- 
tization forces one to choose a special location r™ inside 
the torus which breaks the translational symmetries. 

A possible counter argument is raised: Could TSB be 
avoided by somehow smearing the monopole charge in- 
side the torus? This would presumably restore transla- 
tional symmetry, at least for the direction in which the 



embedding has axial symmetry. It is somewhat surprising 
that the answer is negative, as the particular construction 
below demonstrates. 

Ady Stern [57] has suggested the following construc- 
tion. Consider a large number N of positive Dirac 
monopoles and A^ — 1 negative monopoles placed on the 
middle circle inside the torus, as shown in Fig. 13. One 
arranges both positive and negative monopoles to be at 
equal distances between them such that their mean den- 
sity is uniform. Let us calculate the magnetic fields and 
their effects on bosons on the torus surface, and then take 
N to infinity. 

For large A^, the distribution of monopoles approaches 
a uniform density with total monopole charge of gm- We 
compute numerically the magnetic field created by A^ = 
20 monopoles situated on rings of radius i?, as given by 
Coloumb's law 



B 



2gr3 



(A3) 



The radial magnetic field penetrating the surface of the 
torus in the plane y = 0, becomes increasingly closer 
to a constant, with decreasing oscillatory component as 
shown in Fig. 14. This behavior is precisely analogous to 
the electric field from a ring with uniform charge density. 

Now let us examine the function Wy{;x), evaluated for 
the same monopole configuration as Fig. 15. By Gauss's 
law, it exhibits discontinuous jumps of size 27r (— 27r) at 
each position where the cross section of the torus at x 
cross through a monopole (anti-monopole). Note that 
Wy{x) exhibits A^ positive jumps and A^ — 1 negative 
jumps, corresponding to the number of positive and neg- 
ative monopoles. Therefore, in order for Wy{x) to be pe- 
riodic in X ^ x + Lx, the continuous part of Wy{x) needs 
to compensate for the extra positive jump. Indeed, as 
can be seen by Fig. 15, between the 2A^ — 1 jumps W{x) 
decreases linearly, as demanded by Eq (Al). 

This increasingly discontinuous function does not con- 
verge to a well defined limit function in the large A" limit. 
However, the physically relevant function which effects 
the dynamics of our bosons of charge q on the surface is 
the unimodular phase function 



C(^) 



^iWy{x) 



Coe*«''^''% xe[0,i. 



(A4) 



which is perfectly continuous and periodic on the circle. 
Here we see that C,{xq) = 1 uniquely defines a special 
position Xf) which breaks lattice translational symmetry. 

We note that the values of C,{x) have physical conse- 
quences on the current distribution. In the ground state, 
the loops in the region of C « 1 feel a weak AB fiux, 
and thus a relatively weak persistent current is induced 
in these regions. Similarly, the persistent currents are 
expected to be maximized around loops with C ~ — 1. 

Thus we learn that a static configuration of monopoles 
leads translational symmetry breaking. However, if one 
considers the possibility of an extended quantum wave- 
function of a monopole, the magnetic field will be in a 



15 



quantum superposition. In this case, translational sym- 
metry can be restored in the entangled state of the matter 



field with the electromagnetic field. 



[1] L. Pitaevskii and S. Stringari, "Bose Einstein Condensa- 
tion", Oxford University Press (2003). 
[2] D. Jaksch, P. Zoller, Annals of Physics 315, 52 (2005). 
[3] A. van Oudenaarden and J.E. Mooj, Phys. Rev. Lett 76 

4947 (1996). 
[4] E. Altman and A. Auerbach, Phys. Rev. Lett. 81, 4484 

(1998). 
[5] R. Fazio R and H. van der Zant, Phys. Rep. 355, 235 

(2001). 
[6] M.P.A. Fisher, P.B. Weichman, G. Grinstein and D.S. 

Fisher Phys. Rev. B. 40, 546 (1990). 
[7] Y. J. Uemura et. al., Phys. Rev. Lett. 62, 2317 (1989). 
[8] V.J. Emery S.A. Kivelson, Nature, 374, 434 (1995). 
[9] T. Kostyrko and J. Ranninger, Phys. Rev. B 54, 13105 
(1996). 

[10] A. Paramekanti, M. Randeria, T. V. Ramakrishnan, S. 
S. Mandal, Phys. Rev. B. 62, 6786 (2000); H.J. Kwon, 
A.T. Dorsey P.J. Hhschfeld, Phys. Rev. Lett., 86, 3875 
(2001). 

[11] E. Altman and A. Auerbach, Phys. Rev. B 65, 104508- 
15, (2002). 

[12] M. Franz, A.P. Iyengar, Phys. Rev. Lett., 96, 047007 
(2006); L F. Herbut and M. J. Case, Phys. Rev. B 70, 
094516 (2004). 

[13] A. Mihlin and A. Auerbach, Phys. Rev. B 80, 134521 
(2009). 

[14] A.L. Fetter and J.D. Walecka, "Quantum Theory of 
Many-Particle Systems", (McGraw-HiU, 1971), Ch. 6. 

[15] N. H. Lindner and A. Auerbach, Phys. Rev. B 81, 054512 
(2010). 

[16] V. J. Emery and S. A. Kivelson, Phys. Rev. Lett. 74, 
3253 (1995). 

[17] N. H. Lindner, A. Auerbach and D. P. Arovas, Phys. Rev. 
Lett. 102 070403 (2009). 

[18] We note that in the case of half filled fermion bands, 
particle-hole conjugation is obtained by sending /^ — s- // . 
However, due to the anticommutation of fermion oper- 
ators, Hamiltonians can realize this symmetry only on 
bipartite lattices. In contrast, particle-hole symmetry of 
half filled HCB models occurs on any graph. 

[19] A. Auerbach, "Interacting Electrons and Quantum Mag- 
netism", Springer- Verlag (1994), ch. 13. 

[20] A. Auerbach, "Quantum Magnetism Approaches to 
Strongly Correlated Electrons", Lecture notes com- 
piled by F. Berruto and L. Capriotti, arXiv:cond- 
mat/9801294. 

[21] A. W. Sandvik and C. J. Hamer, Phys. Rev. B. 60, 6588 
(1999). 

[22] K. Bernardet, G. G. Batrouni, J.-L. Meunier, G. Schmid, 
M. Troyer, and A. Dorneich, Phys. Rev. B 65, 104519 
(2002). 

[23] V. L. Berezinskii, JETP 34, 610 (1972); J.M. Kosterlitz 
and D.J. Thouless, J. Phys. C 5, L124 (1972). 

[24] A. Auerbach, D. P. Arovas and S. Ghosh, Phys. Rev. B 
74, 64511, (2006). 

[25] D. P. Arovas and A. Auerbach, Phys. Rev. B78, 094508-5, 
(2008). 

[26] F. D. M. Haldane, Phys. Rev. Lett. 61, 1029 (1988). 



[27] C. Wu H., Chen, J. Hu, and S-C. Zhang, Phys. Rev. A 

69, 043609 (2004). 
[28] F. D. M. Haldane, Y. Wu, Phys. Rev. Lett. 55, 28872890 
[29] F. D. M. Haldane and E. H. Rezayi, Phys. Rev. B 31, 

2529 (1985). 
[30] D. P. Arovas, R. N. Bhatt, F. D. M. Haldane, P. B. Little- 
wood, and R. Rammal, Phys. Rev. Lett. 60, 619 (1988). 
[31] R. Tao and F. D. M. Haldane, Phys. Rev. B 33, 3844 

(1986). 
[32] E. T. Whittaker and G. N. Watson, "A Course of Modern 

Analysis", Cambridge Univ. Press (1999). 
[33] D.P. Arovas, unpublished. 
[34] M. P. A. Fisher and D.-H. Lee, Phys. Rev. B 39, 2756 

(1989). 
[35] Daniel P. Arovas and Jose Freire, Phys. Rev. B 55, 1068 

(1997). 
[36] A. Auerbach and D. P. Arovas, Phys. Rev. B 78, 094508 

(2008). 
[37] C.-C. Joseph Wang, R. A. Duine, and A. H. MacDonald, 

Phys. Rev. A 81, 013609 (2010). 
[38] S. Doniach and B.A. Huberman, PRL 17, 1169 (79). 
[39] D. S. Fisher, Phys. Rev. B 22, 1190 - 1199 (1980). 
[40] W. R. Magro and D. M. Ceperley, Phys. Rev. Lett. 73, 

826 (1994). 
[41] J. E. Avron and R. Seller, Phys. Rev. Lett. 54, 259 (1985) 
[42] G. D. Mahan, "Many-Particle Physics", Plenum, NY 

(1981). 
[43] J. E. Hoffman, E. W. Hudson, K. M. Lang, V. Madha- 

van, H. Eisaki, S. Uchida, J. C. Davis, Science 295, 466 

(2002). 
[44] H.-D. Chen, O. Vafek, A. Yazdani and S.-C. Zhang, Phys. 

Rev. Lett. 93, 187002 (2004) 
[45] L. Balents et. al, Phys. Rev. B 71, 144508 (2005), ibtd 

Phys. Rev. B 71, 144509 (2005). 
[46] We note that a finite superfiuid stiffness has been ob- 
tained for the CB phase by H. Nordborg and G. Blatter, 

Phys. Rev. Lett. 79, 1925 (1997). 
[47] M. V. Feigel'man et. al, Phys. Rev. B 48, 16641 (1993) 
[48] A. S. Sorensen, E. Demler and M. D. Lukin, Phys. Rev. 

Lett. 94, 086803 (2005). 
[49] M. Hafezi et. al, Europhys. Lett. 81, 10005 (2008). 
[50] C. Lannert, Matthew Fisher and T. Senthil, Phys. Rev. 

B 63, 134510. 
[51] Z. Tesanovic, Phys. Rev. Lett. 93, 217004 (2004). 
[52] A. Mehkyan and Z. Tesanovic, Phys. Rev. B 71, 214511 

(2005). 
[53] H. Takagi, T. Ido, S. Ishibashi, M. Uota, and S. Uchida, 

Y. Tokura, Phys. Rev. B 40, 2254 (1989). 
[54] H. Y. Hwang, B. Batlogg, H. Takagi, H. L. Kao, J. Kwo, 

R. J. Cava, J. J. Krajewski, and W. F. Peck, Jr., Phys. 

Rev. Lett. 72, 2636 (1994). 
[55] F. Gerbier, J. Dalibard, "Gauge fields for ultracold atoms 

in optical superlattices" , arXiv:0910.4606vl. 
[56] P.A.M. Dhac, Proc. Roy Soc. (London) A 133, 60 

(1931). 
[57] Ady Stern, private communication. 



